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Abstract 

We study the flux tube thickness of a generic Lattice Gauge Theory 
near the deconfining phase transition. It is well known that the effective 
string model predicts a logarithmic increase of the flux tube thickness 
as a function of the interquark distance for any confining LGT at zero 
temperature. It is perhaps less known that this same model predicts 
a linear increase in the vicinity of the deconfinement transition. We 
present a precise derivation of this result and compare it with a set of 
high precision simulations in the case of the 3d gauge Ising model. 



1 Introduction 



The distinctive feature of the interquark potential in a confining gauge the- 
ory is that the colour flux is confined into a thin flux tube, joining the 
quark-antiquark pair. The quantum fluctuations of this flux tube can be 
well described by an effective string model. The most famous result of this 
model is the well known "Liischer term" which was predicted more than 
25 years ago [1, 2] and was recently observed in high precision montecarlo 
simulations of lattice gauge theories (LGTs) both in (2+1) and in (3+1) 
dimensions with gauge groups ranging from Z2 to SU{N) [3-15]. 

Another well known prediction of the effective string theory is the loga- 
rithmic increase of the width of the fiux tube as a function of the interquark 
distance R. This behaviour was predicted many years ago by Liischer, 
Miinster and Weisz in [16] and subsequently observed in various different 
models, ranging again from SU(2) to Z2 and to the pure gauge percolation 
model [15, 17-22]. Together with the linear rising of the interquark poten- 
tial it has always been considered as one of the distinctive features of the 
confining regime in LGTs. 

A natural question is what happens of the flux tube width as the decon- 
finement temperature is approached from below. According to the above 
picture one would naively expect that the log behaviour should hold in the 
whole confining phase. However we shall show in this paper that this is not 
the case. More precisely, we shall show by means of high precision monte- 
carlo simulations in the (2+1) dimensional Z2 gauge model, that the flux 
tube width also depends on the finite temperature of the theory and that 
near the deconfincmcnt temperature (but still in the confining phase) the 
fiux tube width increases linearly (and not logarithmically!) as a function 
of R. 

As we shall see, this result is in perfect agreement with the effective string 

picture which indeed predicts a logarithmic increase at low temperature but 
shows a much more complex behaviour as the temperature increases and 
ultimately leads, as the deconfinement temperature is approached, to the 
linear behaviour observed in the simulations. 

This paper is organised as follows. In sect. 2 we define the flux tube 
thickness and discuss its evaluation in the framework of the effective string 
model both at zero and at finite temperature. In sect. 3 we present our 
montecarlo simulation while in sect. 4 we discuss our results and compare 
them with the effective string predictions. Sect. 5 is devoted to a few 
concluding remarks. 
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2 Effective string prediction for the flux tube thick- 
ness 



2.1 Definition of tiie flux tube tiiickness 

The lattice operator which is commonly used to evaluate the flux at zero 
temperature through a plaquette p of the lattice is: 

{<t>{p-W)) = ^^-{U,) (1) 

where W denotes a Wilson loop while Up denotes the operator associated 
with the plaquette p. 

In a finite temperature setting we must substitute the Wilson loop with 
a pair of Polyakov loops. The lattice operator becomes in this case: 

/pp't TJ \ 

{<PiP-^P^n) = ^-^^^p^-{U,^ (2) 

where P, P' are two Polyakov loops separated by R lattice spacings. 

Within this setting the flux {(p{p; P, P')) depends on the spacelike co- 
ordinates of the plaquette, on its orientation, on the separation R of the 
Polyakov loops and on the length L of the lattice in the timelike direction. 
It does not depend on the timelike coordinate of the plaquette. Different 
possible orientations of the plaquette p measure different components of the 
flux. In the following we shall neglect this dependence which plays no role if 
one is interested in the R dependence of the flux tube width. Furthermore, 
since we are mainly interested in the tube width half way between the two 
quarks, we restrict the plaquette to lie on the symmetry (hyper)plane half 
way between the two Polyakov loops. Under these conditions we have: 

{<P{p;P,P')) = (^^{h;R,L 

where h denotes the displacement of p from the P P' plane. Along each of 
the directions spanned by /i, the flux density shows a gaussian like shape 
(sec for instance Fig. 2 of [18]). The width of this gaussian is the quantity 
which is usually denoted as "flux tube thickness" : 

w\R,L) = A_ ^ (3) 

Y,j;U{h;R,L)] 
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This quantity only depends on the interquark distance R and on the 
lattice size in the compactified timelike direction L, i.e. on the inverse 
temperature of the model. By tuning L we can thus study the flux tube 
thickness near the deconfinement transition. 

2.2 Effective string model for the interquark potential 

The starting point of the effective string description of the interquark po- 
tential is to model the latter in terms of a string partition function: 



where S^g denotes the effective action for the world sheet spanned by the 
string. In (jH), the functional integration is done over world sheet configura- 
tions which have fixed boundary conditions along the space-like direction, 
and periodic boundary conditions along the compactified, time-like direction 
(the Polyakov lines are the fixed boundary of the string world sheet). 

The simplest and most natural string model is the Nambu-Goto one, 
which assumes the string action S^g to be proportional to the area spanned 
by the string world sheet: 



where g is the determinant of the two-dimensional metric induced on the 
world-sheet by the embedding in and a is the string tension, which 
appears as a parameter of the effective model. 

Eq. (0) is invariant with respect to reparametrization and Weyl trans- 
formations. The standard choice to deal with these symmetries is to choose 
the so called "physical gauge" (see [23] for more details) in which g becomes 
a function of the transverse displacements of the string world-sheet only. 
These displacements (which we shall denote in the following as h^{(;,T) are 
required to satisfy the boundary conditions relevant to the problem — in 
the present case, periodic b.c. in the compactified direction and Dirichlet 
b.c. along the interquark axis direction: 



h\T + L,<i) = h\T,<;)- h'{T,-R/2) = h'{T,R/2) = . (6) 



As it is well known, this gauge choice is anomalouqj: rotational symmetry 
is broken at the quantum level unless the model is defined in the critical 

^Another way to understand this anomaly is to notice that this gauge fixing imphcitly 
assumes that the world-sheet surface is a single- valued function of (r, i.e. overhangs, 
cuts, or disconnected parts are excluded. 






(5) 
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dimension d = 26. However, this anomaly is known to vanish at large 
distances [24], and this suggests to use the "physical gauge" for an IR, 
effective string description also for d / 26. 

If we restrict our attention to the d = 2 + 1 case a further simplification 
occurs since only one transverse degree of freedom (h) survives and the cross- 
interactions terms among different transverse degrees of freedom disappear. 
In the physical gauge, ([5]) takes the form: 

rL/2 fR/2 

S[h] =a dT d?v/l + {drh)^ + {d,h)^ . (7) 

J~L/2 J~R/2 

Since we expect this model to be correct in the large R limit the standard 
way to deal with the square root term in the action is to perform a large R 
expansion in powers of the dimensionless quantity (aRL)^^. 

pL/2 j-R/2 

S[h]r.aLR + - dT dq{Vhf + 0{{aLR)-^) , (8) 

2 J-L/2 J-R/2 

The first term is the classical contribution, it simply gives the area term 
in the interquark potential and we shall neglect it in the following. The 
second term is a standard gaussian action while the higher order 0{{aRL)^^) 
contributions encode string self-interactions. 

If we neglect in the expansion the string self-interaction terms (the so 
called "gaussian approximation" ) , then using standard results of 2d confor- 
mal field theory (CFT) the partition function of the effective string model 
can be evaluated exactly leading to the well known Liischer term. With 
some more effort also higher order terms (and in particular the quartic self- 
interaction term written above) can be evaluated [25] 

The resulting predictions have been compared with Montecarlo simula- 
tions of Polyakov loop correlators for different gauge theories in the last few 
years, showing a very good agreement at large distances and an increasing 
disagreement as smaller distances and/or higher temperatures (i.e. smaller 
values of L) were approached (see for instance [8]) 
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2.3 The effective string width in the gaussian approximation 



The effective string approach aUows to compute the flux tube width accord- 
ing to the following definition: 



where the sum is intended (as in the previous section) over all the surfaces 
bordered by the to Polyakov loops. Since the b.c. are periodic in the timelike 
direction, there is no dependence on t, and on setting x = 0, i.e. half way 
between the two Polyakov loops, we obtain the effective string prediction 
for the flux tube as defined in ([3]). 

The action S[h] should be in principle the whole Nambu-Goto action ([7]), 
but, as anticipated, we shall truncate it to its gaussian approximation ([8]). 
In this way the width becomes a correlator in a 2d free bosonic theory: 



This correlator is singular and must be regularized. The most natural 
choice is a point splitting regularization: 



In fact, in the original lattice description, the flux density is evaluated 
using a plaquette operator which has an intrinsic size of the order of the 
lattice spacing. This ultraviolet scale is translated in the effective string 
description in the e parameter of the point splitting regularisation. 

The correlator in (|2.3p can be evaluated exactly (see the appendix). 
One can then perform an expansion in e of the result. As expected the first 
term diverges logarithmically while the remaining is finite and describes the 
dependence of the result on the modular parameter L/R of the cylinder. 

It is important at this point to distinguish the two regimes: low and high 
temperature. 

• At low temperature, i.e. in the regime in which L » R the flux tube 
width in the gaussian approximation is given by: 




w^{x;R,L) = {h\t,x)h\t,x)) 



w'^{x;R,L) = {h\t,x)h\t + e,x + e)) 




(0)1 



(9) 



q = 
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where (as discussed in the appendix) we use a complex coordinate z to 
describe the cyhnder bordered by the two Polyakov loops (with Re z 
representing the spacehke direction and Im z the timeUke one) which 
are fixed in the positions Rez = Setting Rc = vr|e|/2 we see 

that the dominant term is: 



1 



R 



aw (z) = — log — 



27r 



Rr 



(10) 



as we anticipated, while the next to leading correction in the L » R 
limit turns out to be: 



1 

2^ 



log 



cos 



vr Re z 
R 



which vanishes if we choose z = 0. 

In the opposite regime L « R (i.e. high T, but still in the confining 
phase) the flux tube width has a very different expression: 



aw 



^) = - ^ log ^ + ^ log \0,i2m Re z/L)/e[ (0) | + ^-^^ 
q = e-2-«/i 

(11) 

This can be obtained by a modular transformation of the previous 
result or by direct calculation (see the appendix). In this case in the 
R » L limit the dominant term turns out to be proportional to logL 
instead of log R and the (linear) R dependence only appears in the first 
subleading correction. In fact setting z = and using the relation: 

^l(O) =^2 (0)^3 (0)^4(0) 

we obtain: 

aw\z) = -^iog^-^iog\02mm\ ■ 

Expanding this expression in powers q = e^^'^-^/^ and setting Lc = 
7r|e| we findU: 

O . ^ 1 L R \ n^R , ^ 
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Notice a misprint in the analogous expression reported in [28] where the linear coeffi- 
cient was erroneously quoted to be 1/6 instead of 1/4. 
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This is the major result of this section and we shall devote the next 
sections to a check of this prediction with a set of high precision Montecarlo 
simulations. 

3 Montecarlo simulations 

Testing the logarithmic growth of the flux tube width with Montecarlo sim- 
ulations is a very difficult task since it requires to study very large Wilson 
loops (or Polyakov loop correlators) and to control the statistical errors in- 
duced by the ratio of expectation values of ([2]) 

Both these problems can be solved if one studies abelian LGTs for which 
a duality transformation can be implemented. In particular, for three di- 
mensional LGTs, the dual model turns out to be a spin model. As discussed 
in [18] in this case one can study Polyakov loop correlators of arbitrary size 
by simply frustrating the links (in the dual lattice) orthogonal to the sur- 
face bordered by the Polyakov loops. The expectation value of the energy 
operator (which is dual to the plaquette of the original gauge theory) in this 
environment directly corresponds to the ratio of expectation values of ([2]) 
thus solving at the same time also the second problem mentioned above. 

This strategy was adopted in [18] to study the thickness of flux tubes gen- 
erated by Wilson loops in the 3d gauge Ising model finding a perfect agree- 
ment with the predictions of the gaussian approximation discussed above. 
By choosing different couplings and different Wilson loop sizes the authors 
of [18] were able to test the log growth over a range of more than two orders 
of magnitude. 

The present paper can be considered as a natural continuation of the 
above analysis in the case of the Polyakov loop geometry which, as mentioned 
above, allows to introduce into the analysis also the finite temperature scale 
and allows to study the crossover from a log to a linear growth of the flux 
tube thickness. 

We report here for completeness a few details on the gauge Ising model, 
on the algorithm that we used and on the parameters that we used in our 
simulations. 

3.1 The gauge Ising model 

The 3D TL2 gauge model on a cubic lattice is defined through the partition 
function: 

Zgauge{.ii) = ^ exp (-/JS*) , (13) 
{'Ti=±l} 
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where the action S is a sum over all the plaquettes of the cubic lattice: 



S = -'^an , an = (Ti^cTi^ai^ai^ . (14) 

□ 

This model can be translated into the usual 2>D Ising model by the 
Kramers-Wannier duality transformation: 

ZgaugeiP) OC Zspin{P) (15) 

^ = -^log[tanh(/3)] , (16) 

where Zspin is the partition function of the Ising model on the dual lattice: 
ZspiniP) = Yl ^M-PH{s)) , (17) 

Si=±l 

with: 

H{s) = - Yj{ij)SiSj (18) 
(ij) 

where i and j denote nodes of the dual lattice and the sum is extended to 
the links {ij) connecting the nearest-neighbour sites. For the moment the 
couplings J^jj^ are all chosen equal to +1 . 

Using the duality transformation it is possible to build up a one-to-one 
mapping of physical observables of the gauge system onto the corresponding 
spin quantities. For instance, the vacuum expectation value of Polyakov 
loops correlator can be expressed in terms of spin variables as follows. First, 
choose an arbitrary surface S bounded by the two Polyakov loops; then 
"frustrate" the links of the dual lattice intersecting S, i.e. take J^jj^ = — 1 
whenever (ij) n S 7^ 0. Let us denote with H'{s) the Ising Hamiltonian 
with this choice of couplings: the new Ising partition function Zi.p^j^{/3) = 

^^,^^^cxp (^—f3H'{s)^ describes a vacuum modified by the two Polyakov 
loops, which we shall call the P-vacuum. It is easy to see at this point that, 
thanks to duality we can write the expectation value of the Polyakov loops 
correlator as: 



7' 

z 



/ n exp(-2^s,s,)\ , (19) 



where the product is over all the dual links intersecting S. 
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Similarly it is easy to see that: 



< PP'Wp > 
< PP't > 



(exp(-2/3sfcS/))p 



(20) 



where SkSi is the link dual to the plaquette Up{x) and < >p denotes a 
mean value in the P-vacuum. 

In this way we can immediately obtain the flux density by simply looking 
at the mean value of the (dual of the) plaquette in the model in which all 
the links dual to S are frustrated (see fig [I]). 

More precisely we have: 



3.2 Simulation setting 

We simulated the Ising model (both with and without frustrations) with a 
standard Swendsen-Wang algorithm. We chose to measure the flux with a 
plaquette parallel to the surface joining the two Polyakov loops (i.e. in the 
dual lattice, looking at the product of the spins joined by a link orthogonal 
to such surface, see Fig. We directly evaluated from the simulations the 
flux tube thickness and used a jackknife procedure to estimate the statistical 
errors. We performed all the simulations on a lattice of size 80 x 80 x L. We 
chose f3 = 0.75180 for which the deconfinement transition is known to be 
located exactly at L = 8 [26]. Another reason for this choice is that for this 
value of f3 the string tension is known with very high precision (see [10]) a = 
0.0105255(11). We extracted the flux tube width for values of the interquark 
distance R ranging H from 5 to 50 and for L G {9, 10, 11, 12, 14, 16}, i.e. for 
values of the ratio T/Tc (Tc being the critical temperature) ranging from 
T/Tc = 1/2 to T/Tc = 8/9. 

4 Results 

We report the results of our simulations in Tab. [T] and Tab. [21 The same 
data are plotted in Fig. [2j 

Looking at Fig. O it is easy to see that, in agreement with our effective 
string calculations, the data for L > 9 show a very nice linear behaviour 

^Due to the finite extent of the lattice size in the space directions {L = 80) we expect 
finite size corrections to become important for R > 50. In order to check this expectation 
we simulated for one L value: L = 10 also a few values of R larger than 50. 




(21) 
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Figure 1: Schematic view of the simulation setting. The vertical bonds along the 
horizontal axis represent the frustrated links between the two Polyakov loops. The 
isolated bold face link represents the (dual of) plaquette operator. 
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Figure 2: Flux tube thickness as a function of the interquark distance for various 
values of the inverse temperature L. 
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Table 1: Results of the simulations: the square width as a function of R and L 



L 


R 




L 


R 




L 


R 




12 


5 


20.9(1.6) 


14 


5 


17.9(1.9) 


16 


5 


17.3(1.8) 


12 


9 


34.7(1.3) 


14 


9 


27.7(1.7) 


16 


9 


26.9(1.6) 


12 


13 


44.3(1.3) 


14 


13 


37.8(1.7) 


16 


13 


34.0(1.6) 


12 


17 


59.7(1.2) 


14 


17 


45.6(1.6) 


16 


17 


42.4(1.4) 


12 


21 


69.1(1.2) 


14 


21 


55.0(1.6) 


16 


21 


48.6(1.4) 


12 


25 


83.8(1.2) 


14 


25 


66.7(1.4) 


16 


25 


57.9(1.3) 


12 


29 


95.5(1.6) 


14 


29 


72.4(1.4) 


16 


29 


63.9(1.3) 


12 


33 


107.7(1.6) 


14 


33 


82.2(1.4) 


16 


33 


71.5(1.3) 


12 


37 


122.4(1.6) 


14 


37 


92.6(1.4) 


16 


37 


77.1(1.3) 


12 


41 


131.9(1.6) 


14 


41 


101.6(1.4) 


16 


41 


85.2(1.3) 


12 


45 


144.4(1.6) 


14 


45 


111.3(1.4) 


16 


45 


94.1(1.3) 


12 


49 


160.1(1.6) 


14 


49 


122.2(1.4) 


16 


49 


101.1(1.3) 



Table 2: Same as Tahle\^ 
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as a function of R in the range 10 < -R < 50. The data for ii > 50 show 
deviations due to finite size effects. We performed a hnear fit of the data in 
the range 10 < ii < 50 with the law: 

= k{L)R + c{Ly, k{L) = ^ (22) 

aL 

where, according to the effective string calculation, we should have /cq = 1/4 
The results are reported in Tab. [3l 



L 


k{L) 


c(L) 


^0 


a 




9 


6.19(10) 


23.0(3.2) 


0.587(16) 


4.48(7) X 10"'^ 


2.6 


10 


4.90(4) 


4.1(1.5) 


0.516(5) 


5.10(5) X 10-3 


1.2 


11 


3.85(4) 


2.6(1.2) 


0.446(4) 


5.90(6) X 10-3 
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12 


3.14(4) 


4.4(1.2) 
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6.63(8) X 10-3 


1.1 


14 


2.33(3) 


6.3(1.2) 


0.344(5) 


7.64(11) X 10-3 


0.7 


16 


1.84(3) 


10.6(9) 


0.309(4) 


8.50(12) X 10-3 


0.5 



Table 3: Results of the fit = k{L)R + c{L) for various values of L. In the fifth 
column we also report the values of the string tension a extracted from k{L). 

A few comments are in order on these fits: 

1] We used the following criteria to fix the range of values of R used in the 
fits. In order to fix the upper bound we performed a set of preliminary 
fits for the L = 10 data keeping initially all the data and then itera- 
tively discarding the largest ones looking for an acceptable Xr (i-^- a 
reduced of order unity). In this way we identified as upper bound 
R = 50. For the lower bound we used the same criterion adopted in 
previous works on effective string corrections which assumed the effec- 
tive string to be a reliable description of the interquark potential for 
scales R such that aR? > 1. In our case this means R > 10. Look- 
ing at (jl2p we see that with this choice the first subleading correction 
exp — 27r^ turns out to be negligible within the statistical errors for 
all the values of R and L involved in the fits. 

2] The linear fits show a reduced Xr of order unity for all the values L > 9 
Also the data at L = 9 show a linear like behaviour (see Fig. [2]) which is 
however shadowed by rather large fiuctuations. It is difficult to decide 
if this is the signature of a failure of the effective string picture or if 
the fluctuations are simply due to vicinity of the critical temperature. 
In any case we decided to neglect the L = 9 data in the subsequent 
steps of our study. 
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3] We can compare the values of c{L) extracted from the fits with the 
effective string prediction: 



c(L) = log ^ (23) 



We fitted the data for L > 10 with the law: 

c(L) = alogL + b 



finding a = 17(4) and b = —38(9) with a rather good value of Xr — 
1.6. This result turns out to be in remarkable agreement with the 
effective string prediction: a = l/27r(T = 15.12... (obtained assuming 
a = 0.0105241). 

4] Looking at Tab. [3] we see that the values of k{L) show a L dependence 
different from the one predicted by the effective string. However the 
values that we obtain from the fits smoothly converge toward the pre- 
dicted one as L increase. This can be appreciated looking at the fourth 
column of Tab. [3] where we reported the values of ko extracted from 
the fits assuming a = 0.0105241 which should be compared with the 
effective string prediction /cq = 1/4. For future utility we also reported 
in the fifth column of the table the values which we would predict for 
the string tension if we would fix ko = 1/4 in the fits. These values 
are plotted in Fig. [31 



4.1 Discussion 

The above analysis shows that our data are in substantial agreement with 
the effective string predictions at large values of R and L. At the same 
time however we see an increasing disagreement as R and L decrease. This 
is indeed an usual phenomenon when effective string predictions are com- 
pared with LGT data and points to the well known fact that the effective 
string only represents a large distance effective description of the interquark 
potential. 

We see two main sources of disagreement. The first one, which we already 
discussed above, is the fact that the effective string correction 1/4 in front 
of the k{L) term is only reached asymptotically (see Fig. [3]). The second 
one is that looking at the data in Tab. [1] and Tab. [2] we see that the 
linear behaviour as a function of R extends also in the R < 10 region. This 
definitely disagrees with the effective string prediction which in this region 
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Figure 3: Plot of the string tension a extracted from k{L) as a function of L. 




Figure 4: R dependence of the flux tube width in the R < L region. We also plot 
for comparison a straight line with angular coefficient 1/4. 
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(due to the higher terms in the theta functions) predicts instead a significant 
depletion of the flux tube width (see Fig. (H). 

It is interesting to observe that both these behaviours instead agree with 
a naive Svetitsky-Yaffe [27] dimensional reduction picture (see. [28]). Indeed, 
according to this scheme, the correlation function of two Polyakov loops 
should behave as the spin-spin correlator of the 2d Ising model and the 
combination of plaquette and Polyakov loops used to measure the flux tube 
thickness is mapped into the < eaa > correlator of the same Ising model. 
As discussed in [28] also in this framework one finds a linear increase of 
the (equivalent of the) flux tube width, which however in this dimensional 
reduction scheme holds for all values of R. Moreover the coefficient of this 
linear increase is proportional to 1/m where m is the mass of the equivalent 
2d Ising model and depends on L in a way which strongly resembles the 
behaviour reported in Fig. [3] (see [28]). 

It would be interesting to see if going beyond the gaussian approximation 
one could recover this behaviour also in the effective string framework. This 
would require extending the calculation reported in the appendix to higher 
orders of the Nambu-Goto action or alternatively to adapt to this problem 
the D-brane approach discussed in [29] to deal with the whole Nambu-Goto 
action. Work is in progress to test the feasibility of these approaches. 
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A Appendix A: Derivation of eq.(j9 lllll) 



A.l The Green function for mixed boundary conditions 

The starting point of our analysis is the Green function for a free bosonic 
theory on a cyhndric domain. This Green function is the solution of the 
Laplace equation in a 2d rectangle with periodic boundary conditions in 
one direction and Dirichlet conditions in the other direction. This can be 
mapped to to an electrostatic problem and solved by the method of images. 
We report the result in this appendix for completeness and discuss in some 
detail a few of its properties which will later be relevant for the study of the 
flux tube thickness. 

Let us map the cylinder onto a rectangle of the complex plane centered in 
the origin with sizes [— La;/2, La;/2] x [—Ly/2,Ly/2\. Let us impose periodic 
boundary conditions along the imaginary axis, i.e. for Im {z) = ibLy/2 and 
Dirichlet b.c. along the real axis i.e. for Re (z) = 

Then the Green function can be written as: 

G{z-zo) = -^\og\f{z,zo)\ (A.l) 

with: 

e^[i^{z-zo)/2i^ 

O2 [tt[z + Zq)/ILx\ 
where the Jacobi theta functions 61 and 62 are defined as: 

00 

ei{z) = 2qi J];(-l)"g"("+^) sin(2n + l)z 

n=0 



62{z) = 2q-i ^"^"^^^ cos(2n + 1) 



n=0 

with: 

q = e*^^; r = iLy/2L^ 

In fact, with this definition, log / is analytic everywhere in the rectangle 
except in 2; = zo where f{z) = and the Green function diverges logarith- 
mically. For all the remaining values since log|/| = Re log / (i.e. it is the 
real part of an analytic function) it satisfies the Laplace equation: 

AG{z,zo)=0 
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As for the boundary conditions, by using the transformation properties 
of the theta functions: 



?i(^ + vrT)=-g-V2-0i(z); 



?2(z + vrT)=<7-ie-2-02(2); 



one can immediately see that: 

G{z + iLy, zo) = G{z, zo) 

i.e. that G is periodic along the imaginary axis with period iLy 

In order to check that eq ()A.llA.2p satisfy the Dirichlet conditions along 
the X axis, let us set z = Lx/2 + iy. With this choice we have: 



G{z;zo) = -^log 



[TT{iy - zq)/2Lx + 7r/4] 



O2 [^(iy + zo)/2L^+vr/4] 



from which, using the identity 92{z) = 9i{'k /2 — z) we immediately obtain: 



G{z]zo) 



1 

'2^ 



log 



[^(iy-zo)/2Lx+vr/4] 



Oi [-K{-iy - zq)/2Lx + 7r/4] 



Since the theta functions are real along the real axis, the denominator in the 
above equation is the complex conjugate of the numerator and the argument 
of the logarithm is always unity for any value of y. In a similar way one can 
show that the Dirichlet b.c. hold also for Re 2: = —Lx/2. The expansion 
of eq.s ()A.llA.2p converges very quickly when Ly » L^ ■ In this limit the 
exponentially decreasing terms in G{z, zq) can be neglected and only the 
dominant terms in the theta functions give a contribution. However in the 
opposite limit L^ » Ly such expression is almost useless. In this regime 
one should better perform a modular transformation r ^ — 1/r of the above 
result. In this way one obtains: 



G{z;zo) = -77- log 
zvr 



[iTr{z - zo)/Ly] 



[iTT{z + Zo)/Ly 



+ 



Rez Rezi 







LxLy 



-2nL,/Ly. T- = 2iLx/Ly 



(A.3) 



This expression is equivalent to the above one, but converges well in the 



Lx » Ly limit. 
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A. 2 The flux tube width 

The behaviour of the flux tube width can be extracted from eq. (j2.3|) per- 
forming an expansion in e. Keeping into account the prefactor a in the 
gaussian effective action eq. dSD we see that eq. (|2.3p can be rewritten as: 

aw'^{x,y) = G{z,z + e) (A.4) 

In the hmit e ^ 0, using eq.s ()A.llA.2p and setting = R, Ly = L we 
obtain: 



ayj\z) = _i_ log ^ + i- log (vr Rez/i?) /e[\ 
q = e-V2i? 



(A.5) 



This expression converges well when L >> R. 

In the opposite regime R >> L we must use eq. ()A.3P which leads to: 

aw\z) = - J- log 4^ + J- log |04(27ri Re z/L)/e[ | + 

zvr L zvr LK 

q = e"-2i?/L 
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